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These lectures are meant as an introduction to a few aspects and problems of statistical modeling 
of fully developed turbulence. As such they are clearly not appropriate for experts but they perhaps 
will be more useful for amateurs. 

These series of lectures are essentially divided in two parts. The first part presents quite standard 
materials on fully developed turbulence in Navier-Stokes equation: energy cascade, Kolmogorov 
theory and deviations, etc... The second part deals with the passive scalar problem which aims at 
describing transport phenomena in turbulent flows. This problem provides a pleasant laboratory for 
studying turbulence in which understandable questions may be asked and sometimes also answered. 
It shows intermittency phenomena. The mechanism at the origin of this behavior relies on global 
modes, so-called zero-modes, which are solutions of the inviscid equation of motion of the effective 
field theory. Their existence ensures universality of the intermittency. A similar mechanism may 
hold in Navier-Stokes turbulence. This problem also points towards a few properties of Lagrangian 
trajectories which are at the origin of intermittency phenomena. Namely, the existence of global 
conserved quantities preserved by the Lagrangian trajectories and the breakdown of the Lagrangian 
flow as illustrated by Richardson's law. These properties are expected to be universal in turbulence. 



1 Approaching turbulent flows. 

1.1 Navier-Stokes equation and weak solutions. 

Navier-Stokes and Euler equations, 

Invariance, dimensional analysis and Reynolds number. 

Existence of solutions and weak solutions. 

Since Navier and Stokes fluid motions are suspected to be governed by the Navier-Stokes equa- 
tion which is a differential equation for the velocity field u{x, t) at point x and time t: 

dtu + (n • V) n - vV^u = - (/ - Vp) (1) 

P 

with V the viscosity of the fluid, v ~ 10~^cm^s~^ for the water, p its density, p the pressure and / the 
external force. This equation has to be supplemented by the continuity equation, dtp + V ■ (pu) = 0, 
and the fluid equation of state, f{p,p) = 0. In most cases one considers incompressible fluids for 
which the density p is constant in time and position. We shall set p = 1 by convention. The 
continuity equation then becomes the incompressibility condition: 

V ■u = (2) 

The pressure p is then not an independent variable since choosing the external force to satisfy 
V • / = and taking the divergence of the Navier-Stokes equation (|l|) gives: 

V^p = -(V^n'=) (Vku,) = -(V^Vfe) (n^n,) 

It allows to compute the pressure in terms of the velocity profile. However, the Navier-Stokes 
equation is non-local if it is expressed only in terms of the velocity. 

The Euler equation is the inviscid limit — > of the Navier-Stokes equation in which the 
viscous term —uV'^u has been discarded. It possesses a nice geometrical interpretation as geodesic 
flow on the group of volume preserving diffeomorphisms. 

Once the density has been set to one, the dimensions of the fields are: 

r , length . , length^ length length^ 

M = ^t::— > [p\ = — — 2- ' [/] = T- — 2 ~ 



time time time time 
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The Navier-Stokes equation is invariant under rescalings, 

u{x,t) — > u{x,t) = It~^ u{x/l,t/T) 

p{x,t) — > p{x,t) = l"^ T~'^ p{x/l,t/T) 
f{x,t) ^ f{x,t)=lT-' f{x/l,t/T) 
V — > V = T^^ V 

SO that if u, p, / are solutions of the Navier-Stokes equation with viscosity v so are u, / with 
a viscosity v. It is thus convenient to introduce the dimensionless Reynolds number 

Tie = (3) 

where bi,u is a typical value of the velocity difference on a scale L, a typical length of the system. 
Note that this is a scale dependent definition. For low Reynolds number, He <C 1 the fluid motion 
is regular and laminar. For intermediate Reynolds number TZe of order ~ 1 to ~ 10^ complicated 
flows are observed with, depending on the precise set-up, some of the symmetries permitted by the 
equations of motion and the boundary conditions broken. At higher Reynolds number Tte ^ 1, 
say of order ~ 10^ or more, the fluid motion shows an apparent spatial disorder, which seems to 
be characterized by the proliferation of eddies of all scales, but with statistical restoration of the 
symmetries: statistical isotropy and translation invariance. The limit of infinite, ie. very large, 
Reynolds number is called fully developed turbulence. Formally it corresponds to the inviscid limit 
1/ — > at typical length and typical velocity fixed. 

Fully developed turbulence is carried by irregular solutions of the Navier-Stokes equation. The 
notion of weak solutions was introduced by Leray. It consists in considering solutions in the sense of 
distributions. Namely, u{x, t) is said to be a weak solution of the Navier-Stokes equation whenever 
it satisfies 

J [u^ dt + u^u'' V'^ + uu^ + /^] if^ = , J V = 

for if^ , with V • = 0, and ip smooth functions with compact support. Note that no derivative 
of the velocity field is taken in that definition. This allows to consider weak solutions for which 
the velocity fields possess singularities in its derivatives. Weak solutions to the 3d Navier-Stokes 
equation are known to exist. 

The existence of smooth solutions to the 3d Navier-Stokes or Euler equation with smooth initial 
data is still an unsolved problem. It is known that such solutions exist for short time but their 
existences for all time is still unclear, ie. do they blow up at finite time or not? In 2d the situation 
is different: global existence of solutions to the 2d Navier-Stokes has been proved for a large class 
of initial data. See eg. Q for more informations. Classical textbooks on fluid dynamics and 
turbulence are eg. |Q, ^. 

1.2 Cascades and Kolmogorov theory. 

Inertial range and Richardson cascades. 
Hand-waving Kolmogorov theory and K-ology. 

The Reynolds number @ may be seen as the ratio of the non-linear advection term (u • \/)u 
of the Navier-Stokes equation over the dissipation term —i/'V'^u. So that non-linearity dominates 

2 

in fully developed turbulence. The advection term preserves the energy £ = J ^ whereas the 
dissipation term does not. Prom the Navier-Stokes equation the energy balance is: 

^""^ ^'f-u-HVnf] (4) 



dtj 2 



3 



The first term f ■ u is the energy injected into the system per unit of time while the second 
term z/(Vu)^ is the energy dissipated by unit of time. The pressure does not produce work in 
incompressible fluids. 

In most practical situations the energy is injected into the turbulent system at large scales, 
eg the sizes of the bath in which the fluid moves. It propagates through all system scales and is 
then dissipated at small scales at which the dissipation dominates over the advection. The large 
scale, denoted L, at which the energy is injected, is usually called the integral scale. It serves 
as an infrared cut-off. The small scale, denoted rj, at which dissipation takes place, is called the 
dissipative scale. It serves as an ultraviolet cut-off. The intermediate domain of scales is called the 
inertial range: 

rj <^ Inertial Range <^ L 

The standard characteristic picture for turbulence is a constant transfer of energy from large 
scales to small scales in the inertial range. Following Richardson and Kolmogorov this cascade of 
energy propagates through the scales via a cascade of eddies: big eddies break into smaller ones. 
This picture leads to 'mean field' scaling laws. Suppose that there is a hierarchy of eddies of smaller 
and smaller scale In — ^o?") iT' = 0, 1, • • •, with < 1 the contraction ratio of the eddy size from a 
generation to the following. Let Vi^ be the volume occupied by eddies of size In- Since the density 
of energy carried by the n-th generation of eddies is the total energy accumulated in eddies 

of size In is ^ uf V;„ with Vi^ the volume occupied by eddies of size Z„. The characteristic 
turnover time in the n-th eddy generation is r^^ ~ ln/ui„ and the energy flux at scale / ~ is thus: 

TTl ~ Ei/ti ~ Viuf/l 

Demanding that this energy flux is constant, ie. independent of the scale, vr; ~ const., and assuming 
scale homogeneity in the sense that the volume occupied by eddies of size / is independent of their 
sizes, Vi ~ const., leads to: 

uf el (5) 

with e the energy flux per unit volume, vr; ~ e V/^. Eq.(|5[) is the fundamental scaling law of 
turbulence. It expresses the fact that the energy transfer per unit of time is constant through 
scales in the inertial range. A more precise version of this law, which is known as the Kolmogorov 
4/5-law, will be described in the following section. 
Note that e has dimension 

length^ 

[e] = T — r 

time 

Eq.(^) may be used to estimate orders of magnitude of basic quantities, e.g. the typical turnover 
time, in the inertial range. The rules are that these quantities depend on the scale I only via a 
naive dimensional analysis using eq.(^). Hence, at scale /, the typical velocity is ui ~ gi/3/i/3 g^j^^^ 
the typical time is t; ~ e^^/'^ l"^/^. The typical velocity gradient is 

{Vu)i ~ u/// ~ e^/^r^/^ 

Note that it diverges as / ^ supporting the idea that fully developed turbulence is supported 
by weak solutions, those with possible singularities of the Navier-Stokes equation. This singular 
behavior is only expected in the inertial range, or equivalently only in the inviscid limit of infinite 
Reynold's number. 

The dissipative scale rj may be estimated by looking at the scale at which the advection (u ■ V)u 
becomes of the order of the dissipation —uV'^u. This gives vu^/rf ~ u^/r] and 

„^^3/4--l/4 



4 



As expected it vanishes as v goes to zero, illustrating that r/, or is an ultraviolet cut-off. Below 
this cut-off the advection may be neglected in the Navier-Stokes equation and the velocity field 
becomes more regular. Note that the energy dissipated by unit of time, which may be evaluated 
as — i/(Vn)^ taken at the scale r/, is: 

V (Vn)^ ~ e 

So that the energy dissipation rate is equal to the energy transfer rate, as expected by energy 
conservation. 

The Reynold's number at scale I is TZi = ui l/u so that TZi ~ e^^^ V^l'^ jv. The system Reynold's 
number may be taken as TZi at the large integral scale L, T^e = TLl-, whereas the dissipation scale 
is such that Tty^ — 1. The inertial domain scale with the Reynold's number as: 

As it should be, the inertial range increases with the Reynold's number. 

1.3 Modeling statistical turbulence. 

Random forcing and stationary measure. 
Energy balance and mean dissipation rate. 
More on inertial range and Kolmogorov theory. 

Universality in turbulence is expected to occur in the inertial range only statistically, in the 
sense that statistics of turbulent data extracted from repeated experiments should be more or 
less independent on the precise experimental set-up, eg on the ways the energy is injected in the 
turbulent baths. To theoretically model these repeated experiments one considers random initial 
data and random forcing in the Navier-Stokes equation (^. By universality the inertial range 
statistics should then be independent on the precise statistics chosen for the force. The simplest is 
to choose the force to be gaussian with zero mean and covariance: 

U^{x,t)f\y,s)) = &^{x-y)b{t-s) (6) 

with Vfifix) = to ensure transversality of the force. To mimic the fact that the energy is 
injected at a large scale L, the covariance Cj^{x), which can be chosen in the form C^^{x/L), varies 
on scale L. It is approximately constant up to scale L and decreases exponentially beyond L. One 
may choose C^^{x) to be a gaussian function. 

The velocity field then becomes a random field and one is interested in its multipoint corre- 
lations, and more generally, in its various probability distribution functions. The latter may be 
defined as: 

{5{u{xi,t) -Vi)--- 5{u{Xn,t) - Vn)) 

These are functions of the positions and times at which the velocity is measured. However, one 
expects that at sufficient large time the turbulent systems reach a steady state, independent of 
the initial data. Statistical stationarity means that the equal time correlation functions are time 
independent: 

^(^(Xi,^) • ■■u{Xn,t)) = 

at 

Stationarity implies a balance between the energy injected into the system and the energy 
dissipated by viscous processes. Indeed averaging the energy balance eq.® with = 

yields: 

(/•n) = z.((Vn)2)^e (7) 



5 



The mean injection rate of energy (/ • u) should be equal to the mean dissipation rate i^((Vti)'^). 
As in previous section the energy is injected at large scales of order L and it is expected to be 
dissipated at small scales of order r]: the energy transfer through the scale being without loss. 

The inertial range and energy cascade may then be formulated in a more precise way by intro- 
ducing e<K, the energy dissipated by modes with wave numbers of modulus less than K, and w<k, 
the energy injected into modes of wave numbers less than K. They are defined by: 

e<K = I 7^ I d^xiy{Vu{x)-Vu{{)))e''^ '' 



\k\<K (27r)3 



l\k\<K (2vr)' 

The fact that the energy is injected at large scale means that w<k is approximately constant and 
equal to e as soon as K > 1/L. The fact that the energy is only dissipated at small scale means 
that Ik approximately vanishes at scales K less than l/rj and becomes of order e for K > 1/rj. The 
flux of energy through the surface of wave numbers of modulus K is the difference: 

7f<K = w<K - e<K 

The cascade picture then means that Tf<:K is approximately constant for K in the inertial range, 
1/L <C K <C 1/r]. Looking for constant energy flux is an experimental signal for fully developed 
turbulence. 

The injected energy w<k niay be linked to the Fourier transform of the force covariance (^. 
One has: 

w<K = l;[ 7^4o I d^xe-'^-'^iTCLix) 
2 J\k\<K (27r)'^ J 

So choosing the energy injected spectrum as described above specifies the choice of the force co- 
variance. 



' Energy 
I Injected 



Energy 
Dissipated 



Inertial Range 



K = wave numbers 



Figure 1: Detailed energy balance. 



Velocity correlations in the inertial range <C |x| <C L are extracted by looking at the short 
distance limit in the inviscid correlators. Namely, to satisfy <C |x| one first takes the inviscid 



6 



limit ^ at fixed positions x and fixed integral scale L, and then, one takes the short distance 
limit at fixed L. Assuming the existence of scaling laws for the n-point correlation functions means 
assuming the existence of exponents such that the following limits exist: 

lim lim \~^'' {u{Xxi,t) ■ ■ ■ u{\xn,t)) (8) 

Remark that since in the inviscid limit there is no other scale than L and the positions, taking 
the last short distance limit is equivalent to sending the infrared cut-off to infinity. Inverting the 
limit order would not give the same answer as it would describe correlations in the dissipative 
domain. Since v plays a role of ultraviolet cut-off, velocity correlations behave very differently at 
finite viscosity and in the inviscid limit: they are smooth at short distance at u finite but their 
derivatives develop ultraviolet singularities in the inviscid limit. 

In his seminal 1941 paper [Q], Kolmogorov made a crucial step by postulating the universal 
character of the velocity statistics in the inertial range. In a simplifying way, Kolmogorov scaling 
theory, known as K41, but not to be confused with K-theory, is based on the assumptions that 
there is a constant energy cascade through some inertial range and that the inertial range velocity 
correlators are universal in the sense that they do not depend on how the energy is injected, 
ie. on the force, but only on the positions at which the velocities are taken and on the mean 
energy transfer rate e. As a consequence, the structure functions = ([(n(x) — n(0)) • x]") are 

completely determined by dimensional analysis up to constants: 

Sn{x) ^ {[{U{X) - um ■ XD = Cn {Z \x\r'^ (9) 

Compare this with eq.@. Universality implies that the constants Cn are independent of the forcing. 
The K41 theory also gives a prediction for the energy spectrum E(k), 

E{k) oc e2/3 k-^/^ (10) 

with E{k) the energy in modes of wave numbers of modulus k, 

E{k) dk= [ / d?x {u{x) ■ n(0))e-^'^-^ 

J\k\=K {2TTy J 



Eq.(^) is exact for n = 3, see eq.(p^) below, but the scaling it implies for the higher order 
structure functions seems to be invalidated by real and numerical experiments. 

1.4 Correlation functions and scalings. 

Exact three point function. 

Anomalous scaling, universality and the large L limit. 

The simplest exact result on turbulence is a formula for the order three structure function which 
supports Kolmogorov scaling. At it should be clear from previous sections it is a consequence of 
constant energy transfer. It follows from the Navier-Stokes by imposing the stationarity of the 
velocity two-point function: dt{u{x) ■ u{y)) = 0. This is an exercise worth doing in details. 

Consider the Navier-Stokes equation with random forcing with covariance (P). It implies that 

i-t+St 

u{x,t + 6t) — u{x,t) = / dsu{x,s)ds 

t+st 

2\ 



-, rt+Ot 

u ■ V)u + uV'^u - VpJ (x, t)5t + J fix, s)ds + 0{5t 
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Since the force is delta correlated in time, the term linear in / is of order Jt^/^. Inserting this 
formula into {u{x, t + 5t) ■ u{y, t + 5t)) and keeping term of order 5t only gives: 

dt{u{x, t) ■ u{y, t)) = ([-(n • V)u - Vp] {x, t) • u{y, t)) + i^{V\{x, t) ■ u{y, t)) + {x ^ y) 

+trCL(a; - y) 

The term proportional to {f{x,s) ■ u{y,t)) vanishes since by causality u{y,t) is independent of 
f{x,s) for s > t. The term with the pressure vanishes by incompressibility and others may be 
simplified using translation invariance: 

{Vpix,t)-u{y,t)) = -{p{x,t)V ■u{y,t)) =0 
iy{V^u{x, t) ■ u{y, t)) = -u{Vu{x, t) ■ Vu{y, t)) 

{{u ■ V)u{x, t) ■ u{y, t)) + {x^y) = -^V^ • {{u{x, t) - u{y, t) {{u{x, t) - u{y, t)f) 

Hence stationarity of the equal time two point function gives (we drop the time label to simplify 
the notation): 

- ^V^. • ((u(x) - u{y) {{u{x) - u{y)f) + v{Vu{x) ■ Vu{y)) = ^trC^Cx - y) (11) 

We shall now take the limits v ^ and x — > y in two different orders. Take first the coincident 
point limit at finite viscosity. The first term in eq.(|ll]) then vanishes since velocity correlations are 
smooth at u finite. Sending u to zero afterwards we get: 

lim lim u{Vu{x) ■ Vu{y)) = -trCL(O) = e (12) 



Taking the limit in the reverse order, first i/ ^ and then y ^ x, the second term in eq.(ll) 
vanishes and we get: 

- ^ lini lini V^. • {{u{x) - u{y) ((n(x) - u{y)f) = e (13) 



Assuming isotropy this implies that at short distance 



{{u\x)-u\^)){ui{x)-u\^)){u^{x)-v!'{^))) = -— e {s^ir^ + d^^r' + d^^rA (14) 

15 ^ ^ 

with r = \x\. Equivalently, 

S^{x) = {[{u{x) - n(0)) • xf) = -\-er (15) 

5 

This is known as Kolmogorov four- fifths law. As it is clear from eq.(|l^), e is the mean dissipation 
rate. Remark that the Kolmogorov four-fifths law essentially follows from the cinematics underlying 
the Navier-Stokes equation. 

Of course eq.(|l^) matches Kolmogorov scaling (^) with C3 = —4/5. In dimension d, it would 
be — J^2) • However, there are strong experimental as well as numerical evidences that other 



structure functions do not follow Kolmogorov scaling but 

Sn{x) oc r^" 

with ^„ 7^ n/3. Sample of available data are ^4 ~ 1.28, ^5 ~ 1.53, S^q ~ 1.77, etc. ^2n; 2n > 3, 
should be less than 2n/3. See e.g. for a discussion of the experimental data. 



8 



In the inviscid limit, the only dimensionfull parameter on which the inviscid correlations may 
depend are the positions x, the integral scale length L and the mean dissipation rate e. These 
correlations may also a priori depend on the shape of the force covariance via e~^C^^{x / L) but this 
is a dimensionless parameter. So on dimensional ground, the n-point inviscid correlations are of 
the form 

(e|x|)"/3F„(x/L). 

Kolmogorov hypothesis may then be seen as assuming that the large L — > oo limit exists whereas 
breakdown of Kolmogorov scaling on contrary means that Fn{x/L) diverge for large L. In such 
case the structure functions behave as: 

5„(x) =^„(e|x|)'^/3(|^|/_^)C„W3 as L ^ oo (16) 

Recall that < n/3 so Sn{x) grows with L. Universality is expected to be restored in the sense 
that the anomalous scaling dimensions ^„ are independent of the forcing but the dimensionless 
amplitudes An depend on the details of the forcing via its dimensionless covariance and are therefore 
not universal. Since ^„ < n/3 the higher moments ( |l6|) of the velocity differences become larger 
and larger than their gaussian values as one goes to shorter and shorter scales. This means that 
the probability to have large velocity differences increase at short scales. This phenomena is called 
intermittency. 
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2 Deeper in turbulence flows. 



2.1 Turbulence in 2d. 

Inverse versus direct cascades. 

Enstrophy dissipation but no energy dissipation. 

Kraichnan's scaling, KK-theory and exact 3-point function. 

A special feature which distinguishes two dimensional from three dimensional fluid mechanics is 
the conservation of vorticity moments in the inviscid limit. The vorticity to is defined by w = eijdiUj 
with eij = —eji, e\2 = 1- It is transported by the fluid since it satisfies 

dtuj + {u ■ V)w - vV'^uj = F (17) 

where F = eijdifj with fj the force applied to the Navier-Stokes equation ([l|). As a consequence, 
(?t / cj"" = in absence of force and at zero viscosity. The second moment is called the enstrophy 
so that the enstrophy density is O = ^o;^. The 2d inviscid Navier-Stokes equation admits two 

2 2 

quadratic conserved quantities: the energy / ^ and the enstrophy / 

As first pointed out by Kraichnan in a remarkable paper this opens the possibility for quite 
different scenario for the behavior of turbulent fiows in two and three dimensions, see also [^]. As 
argued by Kraichnan, if energy and enstrophy density are injected at a scale L, with respective rate 
e and rj^ ~ the turbulent system should react such that the energy flows toward the large 

scales and the enstrophy towards the small scales. As this energy flow is quite the opposite to the 
one involved in Kolmogorov's picture for 3d turbulence, one usually refers to the infrared energy 
flow as the inverse cascade and to the ultraviolet enstrophy flow as the direct cascade. 

The fact that energy has to escape to the large scales may be understood by looking at the time 
variation of the energy and enstrophy in absence of forcing, 

9,/y = -V j {Vuf = j 

In three dimensions, the dissipation rate dtu^ /2 is finite in the inviscid limit so that the limit — > 
is accompanied by an increase of the mean square vorticity. This cannot happen in two dimensions 
because the enstrophy is bounded by its initial value since it can only decrease. The total energy 
is then constant in the limit i/ ^ 0, in particular it is not dissipated at small scales. Furtheremore, 
any transfer of energy as to be compensated by a large transfer of energy towards the large scale 
(small momentum k) in order for the enstrophy to decrease, / dk dtE{k) fc^ < 0, and the energy 
to be conserved / dk dtE[k) = 0. At the same time, the enstrophy will be transfered towards the 
small scale (large k) at which it will be dissipated, since the dissipative term —uV'^uj dominates 
over the advection term u ■ 'Vlo at small scale. 

The energy cascade towards large scales is the inverse cascade. It is characterized by the mean 
energy transfer e. Scaling arguments lead to Kolmogorov's spectrum, with E{k) ~ e^/'^ k^^^^ for 
the energy and {5u){r) ~ (er)^/^ for the variation of the velocity on scale r. 

The enstrophy cascade towards small scales is the direct cascade. It is characterized by the 
mean enstrophy transfer 7]^. Since r]^ has dimension 

iVw] = time"^ 

scaling arguments give the so-called Kraichnan's spectrum with 

Eik)^rjTk-' (18) 
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for the energy and 

{6u){r) c::ir]l/^r (19) 

for the velocity variation. In particular, the velocity field is expected to be much smoother in 2d 
turbulence than in 3d turbulence. The scaling behaviors (|l^ ) and (|l^) cannot be true at the same 
time since the scaling ([T8|) for the energy implies logarithmic behavior for the velocity, {6u)^{r) ~ 



rjw^ r'^logr, whereas the scaling (|19| ) for the velocity implies a more regular energy spectrum, 



E{k) ~ k~^~" with a > 0. One sometimes refers to eqs.(|T8,19) to Kolmogorov-Kraichnan's theory, 
but not to be confused with KK-theory. 

As in 3d one may derive exact result for the three-point structure functions under the hypothesis 
that there is no inviscid energy dissipation, 

lim iy{{Vuf) = 



but there is enstrophy dissipation. 



lim iy{{Vaj) ) = r?^ 



As a consequence, not all correlation functions are expected to reach a stationary regime but only 
the galilean invariant ones. 

The derivation of this exact result is similar to the derivation of Kolmogorov's 4/5 law, so we 
shall be more sketchy. The two point velocity correlations {u{x) ■ n(0)) satisfies the same equation 
of motion as in 3d. Hence, in absence of energy dissipation in the inviscid limit, the mean energy 
increases linearly with time according to: 

2 



dt{u'/2) = e (20) 



So, (^) = It up to a constant, and e = ^trCL(O) is indeed the energy injection rate. This is 
simply the obvious statement that in absence of energy dissipation, and/or in absence of friction or 
other processes by which the energy may escape, all energy injected into the system is transfered 
to the fluid. It is expected to be transfered to the mode with the smallest possible momentum, 
the so-called condensate [^. In particular eq.(|20|) shows that in absence of energy dissipation a 
stationary state cannot be reached although structure functions may converge at large time. This 
is one important difference between 2d and 3d turbulence. 

Assuming that galilean invariant correlations, say {{6u)'^{x)) or {uj{x)uj{0)), are stationary, one 
gets two inviscid equations: 

^V'A{6u'){x) {6uf{x)) = 26 -trCi(x) (21) 

and 

-l^'A{Sn'){x) {5uf{x)) = GL{x) 

with Gl = — V'^trCi. One has Gl(0) = 2?]^. These equations are enough to determine the three 
point structure functions. Eqs.(|2l|) slightly differs from its 3d analogue ( [Tl] ) by the extra term 
proportional to e in its r.h.s. which takes into account for the energy increase. 
In the direct enstrophy cascade, x — > 0, this gives 0]: 



{{5u)l) = {{6u). {6u)l) ~ +^ r3 (22) 
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for the transverse and longitudinal correlations, rj^ is equal to the mean enstrophy dissipation rate. 
Thus, as expected the 3-point velocity functions, which only depend on the enstrophy injection rate, 
are smooth and universal in the direct cascade. Eq.(p2|) may be called the "+1/8 law". 
In the inverse energy cascade, x ^ oo, one gets: 

((Hf) = 3((J^)||Mi)=.+|r (23) 

with e the mean energy injection rate. Of course this gives Kolmogorov's +3/2 law for the longi- 
tudinal statistics. 

Although there are necessarely some deviations from Kolmogorov's 'mean field' theory either in 
the energy spectrum or in the velocity statitics in the direct cascade, they seem to be less pronounced 
in 2d than in 3d. A proposition for logarithmic deviations in the direct cascade was presented in 
fiotl . Two dimensional turbulence has recently been observed in remarkable experiments . These 
experimental data seem to indicate almost an absence of, or at least very weak, deviations from 
Kolomogorov-Kraichnan scaling both in the direct and in the inverse cascade. They also show that 
the probability distribution functions of the even order structure functions are close to gaussian 
distributions. More experimental as well as numerical recent results on two dimensional cascades 
are avalaible in refs.[^]. These give support for very tiny, if any, deviations from Kolmogorov and 
Kraichnan energy spectrum in the inverse and direct cascade respectively. But they also indicate 
that correlation functions don't exactly obey gaussian distributions. These experiments were done 
in cells, the bottom of which induce friction on the turbulent fluids. This friction is necessary for 
reaching a stationary regime since there is no dissipative anomaly in two-dimension. The influence 
of friction on the direct cascade has been analysed numerically in |13| and different scenario have 
been discussed in An attempt to use conformal field theory for constructing zero modes in 2d 



turbulence was described in |15]. 



2.2 Dissipation and dissipative anomalies. 

Stationarity equation again and dissipative anomaly. 
Dissipation for weak solutions. 



As it has been apparent in previous section, fully developed turbulence is supported by velocity 
fields which are non-smooth in the inviscid limit. This may be argued either from the scaling 
argument or using the energy balance (|^). At finite viscosity velocity correlations are smooth 
because and the dissipation term in the Navier-Stokes act as ultraviolet regulators. But as — > 
correlations develop singularities. These are not present in the velocity correlations, which are 
still finite in the coincident point limit (because we may measure velocity moments), but they are 
present in correlations of derivatives of the velocity. This a consequence of the fact that velocity 
differences 6u{x) between two neighbour points separated by a distance x scale as 6u{x) ~ x^^^ in 
Kolmogorov approximation. This is also a consequence of the stationarity condition which imposes 
the dissipation to be non zero, 

lim i^iiVuf) = e 

showing that gradients of velocities are diverging at coincident points in the inviscid limit. More 
generally, the dissipation field defined as 

e'^ix) = lim i^{V''u'){x){V^u^){x) (24) 
is non vanishing in the inviscid correlators. This fact is called the dissipative anomaly. 
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By manipulating the Navier-Stokes equation as in previous section, see eqs.(^JT^), one may 
derive an alternative expression for the dissipation field directly as an operator in the inviscid 
theory: 



tre(a;) = — - lim lim V^; • {{u{x) — u{y) {{u{x) — u{y))'^) 

4 x^y u^O 



(25) 



This relation may be thought as some part of the short distance operator product expansion of 
velocity fields in the inviscid limit. Again it illustrates the fact that fully developed turbulence is 
supported by non-smooth alias weak solutions of the Navier-Stokes equation. 

To illustrate the occurrence of the dissipative anomaly let us write the stationarity equation for 
arbitrary correlation functions (On F-niuixn)]) with F„ any functions of the velocity u{xn) without 
derivatives. It reads: 



5uj^{xk) 



+ e' 



k,l 



5ui^{xk)5uji{xi) 



Suii^{xk)Suj^{xk) 

Y[ Fn[u{Xn)]) 



n Fn[u{Xn)]) 



(26) 



The particular case with F[u] = u reproduces eq.(|l3|). The presence of the dissipative field in this 
equation shows that the set of relations satisfied by the velocity correlations as a consequence of 
the stationarity do not form a close set of equations. Eq.(p^) are some kind of equations of motion. 
Note that the precise form of the r.h.s. follows from the hypothesis that the force has a gaussian 
statistics. 

If there are universalities of inertial range velocity correlations, these correlations should not 
depend on the choice of the force statistics, and thus should be characterized by the l.h.s. of the 
equations of motion (26). As we shall see this is the scenario we shall encounter in Kraichnan's 
model: the inertial statistics will be dominated by zero mode solutions of the l.h.s. of the stationarity 
equation. 

Let us do this exercise in details for the two-point functions {Fi[u{xi)] F2[u{x2)]) ■ The gener- 
alization to an arbitrary number of points is simple. It is described in a more abstract way in the 
following section. The stationarity condition for the two-point function is: 

dt{Fi[u{xi)]F2[uix2)]) = 
Using the Navier-Stokes equation and the chain rule, dtF[u] = (dtu) ■ dF[u], it becomes: 

{[-{u ■ V)u^ - V^p] (xi) 6jFi[u{xi)] F2HX2)]) (27) 
+ u{V\^{xi) 5jFi[u{xi)] F[u{x2)]) + {f{xi) 5,Fi[u{xi)] F[u{x2)]) + (xi ^ X2) = 

with 6ji^F[u] the derivative of F with respect to u^''{xk)- The terms in the first line do not pose 
problems and are well defined in the inviscid limit. The second term, which would vanish naively in 
the inviscid limit, is actually non zero due to the dissipative anomaly. Using translation invariance 
to distribute one of the gradient operator on the function F, it may be rewritten in the inviscid 
limit as: 



-u{V\^{xi)6jFi[u{xi)]F[u{x2)) 
i^{{Vu\xi) ■ Vu^{xi))5i5jFi[uixi)]F[u{x2)) 
{e'^{xi)6i6jFi[u{xi)]F[u{x2)) 



uVtv':,{Fi[u{xi)]F2[u{x2)]) 



Xl ^ X2^ 
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with e*-'(x) defined in eq.(24). We used the fact that the second term in the second Hne vanishes in 
the inviscid hmit because the correlation {Fi[u{xi)] F2[u{x2)]) and hence its derivatives are finite 
at = 0. We now have to compute the term in ( p7[ ) involving the force. Since the force is gaussian 
with covariance (pi), functionally integrating by part gives: 



{fix) G[u{y)]) 



SzCfix 



6p{z) ' 



for any functional G of u. So we have to evaluate the functional derivative in the equal time 

limit t = s. By causality this derivative vanishes for s > t. Derivating the Navier-Stokes equation 
shows that it satisfies a first order differential equation 



[dt+M[u]{y)] 



5^^5{y - z)5{t - s) 



5f\z,s) 

for some differential operator M[u]{y) depending on u. As a consequence. 

5u\y,t)] 



e{t-s)G^''{z,y\t,s) 



with 6{t) the step function and G the solution of the linear equation [dt + J\f[u\ (y)] G^^{z, y\t, s) = 
with the initial condition G^^{z, y\t, s)t= 
we get: 



5^^5{z — y). Taking the equal time limit, with 0(0) = ^, 



5u^{y,t)] 



5f\z,s) 



t=s 



y) 



(28) 



This gives, 



{fix) G[u{y)]) = -Cf{x-y) {5kG[u{y)]) 



(29) 



Choosing 0(0) = ^ is a little bit arbitrary but any other choice will give the same result in eq.(27), 
assuming translation invariance. Gathering everything into eq.(27) gives the stationarity equation 



2.3 Fokker-Planck equations. 

Fokker-Planck equation. 
MSR formalism. 



The Navier-Stokes equation with a random forcing belongs to the class of stochastic equations 
of the form 



dtq = V{q) + / (30) 

with V{q) the forces applied to the dynamical variables q and / gaussian noices with covariance, 
{f{t) f{s)) = r 5{t—s), with r some symmetric matrix. For the Navier-Stokes equation the variables 
q are the velocity fields u{x,t). The above method to derive stationarity equations applies as well 
to these stochastic equations. Time derivative of equal time correlation functions are: 

dtiiiij) = E(^(«.-) n i^)+T.(fj n <i^) 
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The last term may be integrated by part since the noices are gaussian: 



if Elk) =^''{^11 1>^) 

Furthermore, since eq.(|30|) is a first order causal differential equation on still has that the variation 
Sq{t) is independent of f{s) for s > t and that it satisfies the linear equation: [dt + 5V[q)]5q = df. 
Thus p/fl \ t=s = 5 ^ij, with the convention that ^(0) = ^ as in eq.(p^). Hence, one gets : 



dt{Uli)=j:i^{q,) I[<lk) + \Y.^^AIl Qk) (31) 

j j k^j i,j k^i,j 

These are generalized Fokker-Planck equations. Similarly, for any function F{q), 



5F{q) I .,5^F{q), 
5qi ' 2 



9,(F(,)) = (y^(,)^) + -P'(^_ 



When the force derives from an action such that V{q) = — F ^ for some action S'(g), a stationary 
measure for eq.(|30|) is provided by the Gibbs measure dq exp[— 5(g)]. 

This does not apply to the turbulence problem because the force does not derive from an action. 
The invariant measure is, of course, much more difficult to find. 

There however exists a path integral formalism to compute correlation functions. It is known 
as the MSR formalism Q. It consists in introducing in the path integral constraints imposing the 
equations of motion. Writing the stochastic equation symbolically as in eq. (|30|) , it leads to the 
following effective action for two sets of variables (p and q: 

S = ij dtip-{dtq-V{q)) + ^J dt{ip-T-ip) (32) 

The variables (p are introduced as Lagrange multipliers imposing the equations of motion. This 
applies to the Navier-Stokes equation with a random gaussian forcing. However it does not provide 
an efficient way to determine the stationary measure, because there is no known efficient way to 
deal with the non- linear advection term which is dominating in the turbulent regime. 



2.4 Multifractal models. 

velocity 6ui{r) ~ uq{1/Iq)^ on domain of dimension D{h). 
anomalous exponents. 

examples: /3 and bifractal models, saturation ^oo- 



Intermittency means that the probability for having large velocity fluctuations increases at short 
scale. Multifractal models are phenomenological models to encode the fact that only part of the 
fluid participates to the cascade: the amount of the participating fluid modes decreasing with the 
scale Although these models do not explain the origin of intermittency phenomena they present 
a simple way to link anomalous scalings to geometrical properties of the fluid behavior. 

Let us first consider the simplest of the multifractal models called the /3-model. In the spirit of 
Richardson cascade, let us suppose that there is a hierarchy of eddies of scales In = ^o?") = 0,1,- ■ ■, 
with ^ < 1. Now suppose that, contrary to the usual Richardson cascade, the volume V;,^ occupied 
by the eddies decrease from one generation to the next by a factor (3 so that 
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with 3 — D = log /3/ log? > 0. Demanding a constant energy flux e = m "^Viu^ /I yields 

5ui ^ no (///o)^''"'^/' 

for the typical variation velocity 5ui at scale /, with Uq = IIq. Assuming that the probability to 
find eddies of size / at some point in the fluid is proportional to the volume occupied by them, the 
moments of the velocity variation may be estimated as 

{i6uir) ~ ul {l/h)^- with = p/3 + (3 - D){l-p/3) 

which describe some deviations from Kolmogorov's scaling. 

More generally, assume that for x in a subset Vh with dimV/i = D/^ the velocity variation 
5ui{x) = u{x + 1) — u{x) scales for / ^ as 

5ui{x) uq for X G Vft, d\mVh = Dh (33) 

The possible scaling dimensions h are supposed to belong to some interval /imin < h < /imaxi and 
Dh form the spectrum of fractal dimensions. 

To compute the moments of the velocity variation one then assumes that the probability to 
find the scaling behavior ( p3| ) for the velocity variation 6ui{x) is proportional to probability for the 
segment between x and x + / to cross a point in Vh- This probability is proportional to the volume 
of Vh thickened along its transverse directions on a depth of order /. 




Figure 2 : Thickened fractal volume. 



This volume scales as ~ l'^ ^'^ with d the dimension of the ambient space, usually d = 3. Hence, 

with dfi{h) some measure encoding the probability distribution of the fractal volumes Vh- In the 
short distance limit, / — > 0, the above integral is dominated by a saddle point and thus 

{{Suif) ~ < (^^) with = min {ph + d- Dh) (34) 

Since intermittency means that velocity variations become larger at short scales, one expects that 
the spectrum of h is such that h < /imax = 1/3. Of course < Dh < d so that < p/3 as expected. 

The relation between the fractal dimensions Dh and the exponents is a Legendre transfor- 
mation. It can be thus inverted as 

Dh = min {ph + d - ^p) (35) 
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As it is clear, this geometrical construction does not predict the spectrum of anomalous dimensions. 

The /3-model is a fractal model with a single fractal dimension h = {D — 2) /3. The next simplest 
multifractal models are bifractal models with h taking only two values /iq and hi. This is a case we 
shall encounter in one dimensional Burgers turbulence in which ho = with D{0) = (ie. points) 
and hi = 1 with D(l) = 1 (ie. intervals). The anomalous dimensions are then = min(l,p). It 
saturates to the value = 1 as soon p is bigger than one. 

The fact that the anomalous dimensions saturate to some finite value ^oo when p — > cxo as 
an interesting geometrical interpretation. Indeed, assuming that h belongs to some interval so that 
Dh is bounded, then finiteness of ^oo is compatible with eq.(]3^) only if /imin = 0. From eq.(]33|), this 
means that there exists a subset of the fluid volume Vq of dimension Dq = d — ^oo on which the 
velocity is discontinuous, since 5u{x) ~ 0{1) for x G Vq. In other words, finiteness of ^oo implies 
velocity shocks. 

Multifractal models may also be used to estimate behaviors of other quantities such as the dissi- 
pation. The inviscid relation (25) between the velocity and the dissipation field, e(x) oc di{5iu{x))^ , 



tells us that for x on a subset Vh of dimension the dissipation scales as: 

6ie{x)^e{l/lof''~^ for x G Vh 

Here, Si e(x) may be defined by averaging e{x) on a ball of radius / around x. As it should 5i e(x) 
is singular as / ^ 0, for /imax = 1/3- Their moments behave as 

{{dief) c^ePil/kr" with 7]p = mm{{3h-l)p + d-Dh)=C3p-p (36) 

h 

Contrary to velocity correlations, dissipation correlations decrease with the distances in the inertial 
range since < p. 
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3 Scalar turbulence. 



3.1 Transport and Lagrangian trajectories. 

Transport equation. 
Lagrangian flows. 

Definition of the pdf's (backward/forward). 



Scalar turbulence, which aims at describing transport phenomena in turbulent flows |16|, pro- 
vides a simple toy model for modeling turbulence statistically. The equation governing a passive 
scalar advected by a turbulent flow and subject to a small dissipation is: 

dtT + V ■ {u T) - kAT = f. (37) 

Here T(x, t) represents the passive scalar, e.g. the density of another passive fluid. The forcing term 
/ is here to compensate the energy dissipation caused by the term proportional to the molecular 
diffusivity k. The velocity field u with V • u = is supposed to be random. 

The transport equation (|37|) is linear and can thus be solved explicitly. Suppose for simplicity 
that at initial time, T{x,to) = Tq{x), the solution to eq.(|37D is then: 



T{x,t)= I dsdyR^,{x,t\y,s)f{y,s)+ [ dy Rf,{x,t\y,to)To{y) (38) 
Jto J 

with i?K the elementary solution of (5t + V • u — kV'^)Ri^ = with initial condition, Rk{x, toly, to) = 
6{x - y). 

As a consequence of eq.(38), correlation functions of the scalar may formally be written as: 

(H T{xj,t)) = / dsdy (H RniVj, Sj)) (H /(%•, s,)) (39) 

i=l 3 3 

in the case with zero initial condition Tq = 0. The force and trajectory statistics have been 
f actor ized. 

In the limit of zero molecular diffusivity, k — > 0, properties of the passive scalar are of course 
intimately related to behaviors of passive particles advected by the fluid. Consider a passive particle 
initially at position xq at time to in 8' turbulent fluid and transported by it. It follows a trajectory, 
called a Lagrangian trajectory, whose equation of motion is the first order differential equation: 

x{t) = u{x{t),t) with x{to) = Xq (40) 

with u{x, t) the velocity field. When needed to specify the initial data, we shall denote by x(t|xo, to) 
the solution of eq.(pC|). 

In the limit k — > 0, the resolvent i?^ may formally be expressed in terms of solutions of La- 
grangian trajectories: 

RK=o{x,t\y,s) = 6{x - x{t\y,s)) (41) 

with x{t\y,s) the time t position of the particle starting at y at time s. The solution (^) then 
codes the fact at k = the scalar is simply transported by the fluid. It provides an explicitly link 
between properties of the scalar and those of the Lagrangian trajectories. 

Lagrangian trajectories may be described more precisely by introducing their probability dis- 
tributions. One may specify either the initial or the final positions. In the former case, one looks 
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for the distribution of the final positions. This is encoded in the forward probabihty distributions 
which are formahy defined for n particles as: 

P{oTwa.rd(^^ Axo,to) dx = Proba {x{t) G [x,x + dx]; x{to) = xq, to < t) 

n 

= (Y\5{xj - x{tj\xQj,tQj))) dx (42) 



where x(i|xo,to) is the solution of the Lagrange equation (|4(]| ) with initial condition x(to) = 2:0. 
This is the time t position of the particle which was at xq at time to- The average is taken over the 
velocity field realizations. 

Conversely one may specify the final positions and look for the distribution of the initial posi- 
tions. These are formally defined as: 

-fbaikward(^' ^1^0, ^o) ^3^0 = Proba (x(to) G [xq, 3^0 + dxQ\; x{t) =X,t> to) 

n 

= {Y[S{xoj - x{toj\xj,tj))) dxo (43) 
i=i 

with x{to\x,t) the position of the trajectory at time to which will be at x at later time t > to- 

Caution has to be taken when the velocity if not regular enough since, as we shall discuss 
later, in such cases Lagrangian trajectories may not be uniquely defined. This is expected to 
be the case in the inviscid limit i/ — > but not at finite viscosity. So in turbulent flows, these 
probability distributions have to be understood as the limit of their regularized analogue obtained 
by considering finite viscosity. For incompressible and time-reversal invariant velocity fields these 
probability distribution functions coincide. 

The velocity probability distribution may be reconstructed from the trajectory probability 
distribution functions. For example, by expanding the trajectory equation (^) to lowest order, 
x{t + e) = x{t) + eu{x{t),t) + • • •, one gets: 

Froha {u{y, t) e [v,v + dv]) dv = {5{v - u{y,t))) dv (44) 

= lim Pbackward(y, t\y - ev,t-e)d (ev) 

This formula applies only if the velocity is finite at each instant and it assumes that trajectories 
are well-defined such that one may expand eq.(Eo|) to lowest order. 



3.2 Kraichnan's passive scalar model. 

Definition, motivations and properties of the velocity fields, 
Two-point function, energy cascade and inertial range. 



In Kraichnan's model of passive advection, |17], the statistics of the velocity field in eq.( 
independent of the forcing, is supposed to be gaussian with zero mean and with the two-point 
functions 

(n"(x, t)u\y, t')) = D''\x - y)5{t - t') with V aD"^ = 0. (45) 

To analyze scaling properties we shall use the following expression for D""^: D"'^{x) = D{Q)5°'^ — 
d^-^ix) with 



d''\x) = D[{d + i-l)6''^ -i'^]\x\^ (46) 
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where ^ is a parameter, < ^ < 2. A more rigorous approach requires regularizing by 
introducing a infrared cut-off. For example: 

r „ik-x ( h"-h^\ 



D''\x) =Dq j dk 



(A;2 + rn2)(rf+C)/2 I A;2 



Clearly, this distribution for u is far from realistic. It mimics however the growth of the correlations 
of velocity differences with separation distance, typical for turbulent flows since 6u{x) ~ r^^'^ . The 
fact that the two-point functions (^5|) are white noise in time is crucial for the solvability of the 
model but it is of course very far from reality. It in particular implies that, as in Brownian motion, 
at each instant all the velocity moments t)^) are infinite. 

As for modeling turbulence, the forcing term is also assumed to be gaussian with mean zero 
and two-point function 

{f{x,t)f{y,t')) = C{^)5{t-t') (47) 

The rotation-invariant function Cl{x) = C{x/L), which could be chosen to be a gaussian, varies 
on scale L. 

The parameter ^ fixes the naive dimensions under rescalings x fix, L fj,L. These are: 
M=C/2 ; [T] = (2-0/2 ; M = (2-e) 

The value ^ = 4/3 corresponds to Kolmogorov scaling. Indeed when comparing with Kolmogorov 
scaling one has to remember that Kolmogorov typical turnover is finite and scales as — r^/^ so 
that (((5n)2(r)) ~ r^/^ ~ r^/^/r^. This has to be compared with the two point function {{5u)'^{r)) ~ 
6{t) in Kraichnan's model. 
Kraichnan's model provides a toy model for turbulence with a scale domain, the inertial range, 
with constant energy flux. The scalar energy density is £ = ^T^. In presence of non-zero molecular 
diffusivity, the mean energy balance reads: 

dtj£ = j [{f-T)-K{{VTf) 

So, as in previous section, existence of a stationary regime requires 

6 ^K((Vr)2) = (/.T) = 1(7^(0) (48) 

with e the mean dissipation rate. As in developped turbulence, one has a dissipative anomaly since 
the mean dissipation rate k((VT)2) does not vanish in the limit k — > 0. Eq.(|48|) simply means that 
the amount of energy transfered to the scalar by the force equals the amount of energy dissipated. 

Stationarity equation may be derived as in previous section, see eqs.( |Tl|j26|) . Again let us do it 
in some details. For example one may expand T{x,t + 5t) to leading order and insert it into the 
correlations functions to derive the stationarity equations. Using eq. (|38|) one has to expand the 
resolvant which is formally define by the path ordered exponential: 



Ri^{x,t\y,to) 



P-exp[-/ ds (kV^ - V • u)(s) 

Jtn 



(49) 



To leading order in 6t this reads: 

t+5t t+St si 



R^{-,t + 6t\-,t) = 1+ j dsV-u{s)-KV^6t+ j dsi j ds2{V ■u){si)(y ■u){s2) + 
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Recall that since both / and u are delta-correlated in time the first term in the r.h.s. is of order 
but the last term is of order {5t). As a consequence one has for the two point functions: 



(r(x, t + 5t)T{y, t + 6t)) - {T{x, t)T{y, t)) = k (V^. + j {T{x, t)T{y, t)) 6t 

t+St si 

+ J J '^'^2 ((V • n)(x, si)(V • u){x, S2)T{x, t)T{y, t)) + (x ^ y) 
t t 

t+St t+St 
+ J dsids2{{V ■u){x,si){V ■u){y,S2)T{x,t)T{y,t)) + J dsids2{f{x,si)f{y,S2)) 
t t 

All correlation functions may be factorized since by causality T{x,t) is independent of u{y,s) or 
f{y,s) for s>t; eg. {u{x, si)u{y, S2)T{x,t)T{y,t)) = {u{x, si)u{y, S2)){T{x,t)T{y,t)) for si,S2 > t. 
Assuming translation invariance, stationarity of F2{x) = {T{x)T{0)) then gives: 

(_^V2 - d'''{x)VaVb) F2{x) = Cl{x) (50) 

The fact that it is a close equation reflects the linearity of the transport equation and the particular 
choice of the velocity statistics. This property makes Kraichnan's model theoretically attractive. 
In the rotational invariant sector the differential operator = — d"^(a;)V"V^ is 

Ml = -D{d-l)^^r''+^-'4- 
With appropriate boundary conditions, the solution to eq.(^0|) is: 

^^(^) = X Did - + 
It possesses an inviscid limit k ^ whose short distance expansion is: 

The above const., which may be read of in eq.(^), depends on the details of the shape of the 
forcing correlation Cl{x). It is thus non universal. It however cancels in the two-point structure 
functions 

5,(.) = ((T(.) - mf) = r'-' + ■ ■ ■ (53) 

It is universal since it only depends on the distance and on the mean dissipation rate. Note that 
we recover the naive scaling dimension of the scalar [T] = (2 — ^)/2. 

To pursue further the analogy with the turbulent cascade, we should exhibit a constant energy 
transfer in some scale domain. As in previous section, we introduce the energy dissipated £<k into 
modes of momenta k less than K, as well as the energy J-<k injected into these modes: 

£<K = f [d'xK{VT{x)VTme-^'-^ 

j\k\<K {^T^r J 

^<K = I 7^ f d\{T{x)fme-^'-^ 

J\k\<K [ZTT)'^ J 
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Fourier transforming eq.(^) then gives the energy balance 

Vkk = J^<K — £<K 

with V<K the energy flux through modes of momenta k, \k\ = K. Using the expression ( |5T|) for 
one then verifies that these quantities have the following K dependence in the limit of zero 
molecular diffusivity, similar to the real turbulent case: 



Energy 
Injected 



Energy ■ 
Dissipated / 



Inertial Range K = wave numbers 



Figure 3: Scalar energy balance. 



Hence the energy transfer V<k is approximately constant in the domain rj ^ x <^ L with 
7/ ~ {k/D)^/^ the dissipative scale below which dissipation dominates over advection. This is 
similar to fully developed turbulence. 

3.3 Anomalous scalings and universality. 

Two-point function with friction. 
A'^-point functions and zero modes. 
Anomalous scaling and universality. 

If the scalar correlation functions were governed by a mean field theory, the structure functions 
would scale as ((r(x) - r(0))^) ~ (er)^(2-0/2_ instead, we shall see that in the inertial range the 
even structure functions have anomalous scaling with [18|: 

/ r \ Pat 

((r(x,t)-r(o,t))^) ^ A^l^^j l^|(2-OA^/2 + ... (54) 

The anomalous exponents have been computed in a ^-expansion or l/d-expansion ||l^, |l9|: 

The exponents are universal depending only on ^ but the amplitudes Aj^ are not: they depend 
on the shape of the covariance Cl. These anomalous exponents are the signal of intermittency 
phenomenon, with large fluctuations at short distances. 
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Let us now return to the problem of computing the A^-point correlation functions in Kraichnan's 
model. First we have to find the constraints imposed by the stationarity condition 



dt{T{xi)---T{xN)) =0. 



Using either ones of the methods explained in previous section leads to the following stationarity 
equations: 



1 



(r(xi) • • • T{xn)) = -J2 ^Lixjk) (T(xi) • • • T{xj) ■ ■ ■ T{xk) ■ ■ ■ T{xn)) 
where overhatted quantities have to be omitted and, assuming translation invariance, 



(56) 



(57) 



Equations ( p6| ) form a triangular set of equations which recursively determine the A^-point functions. 
In the inertial range, scalar correlation functions satisfy eq.(p^) but with k — > 0, ie. with 



M% = Wd^\x,,)Vt^Vi^ 



(58) 



This is a homogeneous singular second order differential operator which possesses homogeneous 
zero modes that we shall denote (pn{x): 



M% ipn{x) = , (/J„(Ax) = ipn[x) 



(59) 



No explicit expressions for the zero modes are known, only ^ or \/d expansion have been found, 
see eq.(^. 

As in the case of the two-point function with friction described below, one can then show that 
these zero modes provide the dominating contributions to the inertial range correlation functions: 



(r(xi) • • • T{xn)) = An Mx) + 



(60) 



with pjv = (2 — 0-^/2 — ^N- The dots refer to subleading terms as L — > oo. The anomalous scaling 
behaviors announced in eq.(|5^) then arise from these zero mode contributions. 

It is important to realize that universality of the scaling laws (^4|) is ensured by the fact that 
A'^-point functions are dominated by zero modes. The scaling dimensions only depends on the 
velocity field but not on the force statistics, even if the latter would not be gaussian. 



The large A^ limit of the anomalous dimensions ^at have been estimated numerically [21| or 



using instanton computations |22|. One has: 



Saturation value 


X Anomalous 




/ Exponants 




N 



Figure 4 : Scalar anomalous dimensions. 
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So saturates as A'" ^ oo: ^jv — ^oo- Instanton calculations [|2^ give ^oo — si 
d(2 — ^) ^ 1. In a multifractal description, see Section 2.2, this means that the minimum of the 
spectrum of fractal dimensions for the variation of the scalar vanishes: /imin = 0. Thus, there are 
subsets Vo of the ambient space, with dimVo = Dq, on which the scalar varies discontinuously. For 
typical velocity realization, the ambient space is decomposed into domains, on which the scalar 
varies continuously, whose boundaries Vo are shock locations. 

To illustrate the origin of these zero modes and anomalous scalings, let us consider a simpler toy 
problem. It describes a passive scalar advected by a smooth random velocity field but in presence 
of friction ||2^. The friction is represented by an extra linear term, so that the scalar equation of 
motion is: 



dtT + V- (n T) + T/r - kAT = f. (61) 
with T the relaxation time induced by the friction. The velocity is supposed to be gaussian with 



two-point function (45,^) with ^ = 2 such the velocity correlations are smooth functions of the 



positions. At k = the stationarity equation for the two point function F2{x) = {T{x)T{0)) reads: 



T 



This is solved by factorizing F2{r) = H(r)G{r) with G solution of this differential equation without 
r.h.s.. This zero mode solution are homogeneous, G{r) oc with 



d (P 2 

"* = -2^\/T + D(d3T)7 («^' 

With appropriate boundary conditions at infinity and at the origin, the result for F2 is then: 



dp 



p 



At short distance, r — > 0, it behaves as: 

F2(r) ~ const. r'*+ +••• (63) 

Naive scaling analysis would have given scaling dimension zero for T, since ^ = 2. Instead the two- 
point function is dominated by a zero mode solution of stationarity equation and as a consequence 
it acquires an anomalous scaling behavior. Furthermore, the scaling exponent a+ is only function 
of the velocity statistics and it is thus universal, whereas the amplitude depends on the details of 
the force correlation and it is thus not universal. 

More informations on the scalar turbulence problem may be found in [^] and references therein. 
I apologize for all possible omissions of this list. 
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4 Lagrangian trajectories. 



4.1 Richardson's law. 

Explosive separation of trajectories chaos) 
Breakdown of Lagrangian flows. 



Richardson's law describes phenomenologically the separation of passive particles transported 



by turbulent flows [^]. Recall the definition (40) of Lagrangian trajectories: 



x{t) = u{x{t),t) 

with u{x, t) the velocity field. Consider now two particles at positions xi and X2 and let p{t) = 
xi{t) — X2{t) be their relative position. Assuming Kolmogorov's scaling, it satisfies in the inviscid 
limit, 

p{t) ~ e^/^ (g4) 
Thus, p^/^ — Pq^^ ~ e (t — to)- If the particles are at identical positions at to = , then 

p'^it) ~ et^ (65) 

This is Richardson's law: trajectories separate in time polynomialy. It describes the explosive 
separation of Lagrangian trajectories in the limit of infinite Reynold's number, ie. infinitely close 
trajectories separate in a finite time. 

This behavior is different from properties of chaotic systems in which the separation of initially 
close trajectories is exponential. Indeed, suppose the trajectory equation would have taken the 
form X = v{x,t) with v{x,t) smooth enough such that its variation on a scale p is 5v ^ a p. Then 
the analogue of eq.(|^) would have been p a p whose solution is exponential: p{t) ~ poe'^*. In 
particular, if trajectories start initially at identical position, pQ = 0, then they coincide at later 
time, ie. p{t) = 0. It is also clear from this tiny evaluation that the explosive behavior ( |65| ) is 
linked to the non-smoothness of the velocity field. 

Eq.(|65|) means that two particles starting at positions which are indistinguishable in the inviscid 
limit, ie. two initial positions distant less than the viscous scale r/, don't necessarily follow the same 
trajectory. This statement has to be understood in a probabilistic sense since eq.(|65|) is a mean 
field estimate for the averaged separation. Richardson's law indicates that Lagrangian trajectories 
are ill-defined in the inviscid limit since there may be different solutions to eq.(|4^) with identical 
initial data. This explosive behavior may be rooted to the non-smoothness of the velocity field in 
the inviscid limit, since if it were smooth enough then Lagrangian trajectories (^), which are first 
order differential equations, would be uniquely defined by their initial conditions. 



4.2 Lagrangian flows in Kraichnan's model. 

Path integral representation of the pdf's. 
Heat kernel representation. 

Recall equation ( |39|) which links properties of the scalar to those of the Lagrangian trajectories. 
To take into account the diffusion of the passive scalar, we add a white-noise to the Lagrange 
equation and look for the statistics of trajectories defined by: 

x{t) =u{x{t),t) +'q{t) (66) 
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with u{x, t) a gaussian velocity field with statistics ( [dq) and r] a gaussian white-noise with two-point 
function, 

{7l{t)r,{s)) =2K6{t-s) 

This modification will also take into account possible ill-definedness of the probability distribution 
functions of trajectories due to the non-smooth nature of the velocity field in Kraichnan's model. 

Since u{x, t) is time-reversal invariant, the backward and forward probability distribution func- 
tions coincide. They admit a simple representation as a heat kernel for the differential operator 



Pj^'^\x,t\xo,to) = [exp(-|t-to|A^^) 



{x\xo) 



(67) 



with Ai'^ defined in eg. (|57|) . This representation of P^'^^x , t\y , s) together with eq.(|3^) provides 
another way to derive the stationarity equation (|5^). 

A possible derivation of eq.(|67|) uses path integral manipulation. By introducing a path integral 
representation of the delta-function coding the trajectories, one has: 



N 

PK^\x,t\xo,to) = {YlS{Xj - x{tj\xoj,toj))) 

i=i 

The average is over u and r/. Since they are both gaussian we get: 
Pi^) (x,t|xo, to) = 

rt 
'to 



(68) 



(69) 



x(t)=x 



J DpDx exp -Jdsl Y.{kp]{s) + ipj{s) • Xj{s)) + ^ ^ D'-^xj - Xk)p'-{s)pl{s) 

r.\=Xr. . *° \ i 



x{to)=Xo 

This is the path integral representation of the heat kernel of A4^. 

In absence of white-noise in Lagrange equation, i.e. in the limit of k ^ 0, we get: 



Pi=l{x,t\xo,to) = [exp (^-\t - to\M%) 



{x\xo) 



(70) 



with M-% the k = hmit of M^j defined in eq.(p^). For ^ = 0, M% reduces to the Laplacian 
operator and Pj—l to the usual heat kernel of the Brownian motion. 
The kernel ^^=0 satisfy the semi-group law: 

J dyPS(x,0|y,.)PS(2/,.|z,t) = Pi3(x,0|z,t) 

They thus define a Markov process. Translation invariance ensures that p(^) only depends on the 
time difference. Time-reversal invariance of the velocity statistics implies that P^^\x,t\y,s) = 
P^^\y, t\x, s). Due to the homogeneity property of A4% we have 



A^^pS(Ax,A(2-^)t|Axo, A(2"5)to) = pS(x, t|xo, to) 
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4.3 Slow modes. 



Hyperdiffusivity. 

zero modes as conserved quantities and slow modes. 

The zero modes, which dominate the scalar correlations in the inertial range, have a simple 
interpretation in terms of Lagrangian trajectories: They are quantities conserved in average. 

One may test the distribution of the trajectories by computing averages of some functions /(x) 
of the particle positions: 

(/){t[xo,to) = / dxPl^l{x,t\xo,tQ) f{x) (71) 

This is a function of the initial positions and of the time t at which we evaluate it. To avoid 
irrelevant technical problems we shall only consider translation invariant test functions. 

Generically, these averages grow with time indicating an increase of the distance between the 
particles. More precisely, if the test function is homogeneous of degree a, /(Ax) = \" f{x), then 

(/)W.0A.) ^ t''"^^-^^ (72) 

When ^ = one recover the usual diffusive behavior of the Brownian motion. For ^ < 2 the increase 
is faster and one usual refers to it as super-diffusive. Note that the exponents diverge as ^ — > 2. 
Eq.([72|) follows from homogeneity property of the heat kernel -P^=o- 

The zero modes are characterized by completely atypical behavior. If (/?n(x) is a zero mode, 
M^ifn = 0, then: 

(v5n){t|a:o,io) = V'nl^^o), independent of time (73) 

So the zero modes reflect coherent structures which are preserved by the flow. The derivation of 
eq.([73|) is simple. By definition ([7l|), one has: 



9t(9?n)(i|xo,to) = J dxdtPj^^l{x,t\xo,to)fn{x) 

= -{{M%ipn)){t\xo,to)=0 

where we used that {dt + M%)Pj;.^l{x, t\y, s) = 0. 

Besides the conserved zero modes 9?n(x), one actually has a series of modes, that we denote as 
Vn\k{x), k > 0, whose correlation functions increase slower with time than the hyper diffusive law 
They are related by: 

M% ipn-k+l{x) = 9?n;fc(x) (74) 

They form a tower of descendants whose top is the zero modes = ^Pn-fi- In particular, if the 
zero mode has dimension then its descendants ipn;k has dimension ^n;k = + (2 — These 
modes code for the behavior of the Lagrangian trajectories as points become closer. More precisely 
the probability distribution functions of the trajectories admit a short distance expansion of the 
following form: 

pS(x,i|A2/,0) ^Y.^^"''''M^^t)fn;k{y) (75) 

n;k 

as A ^ 0. The descent equation (|7^ follows from consistency conditions for this expansion. 
Indeed inserting the defining relation [dt + A4^)P^^q(x, 0) = into this expansion implies 
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{dt+M%)'ilJn;k{^,t) = 0. Similarly, plugging the expansion (|7|) into {dtPl^^l+Pj^^lM%){x,t\y,0) = 
yields: 

where we used the fact that M.% is homogeneous with scaling dimension — 2). Comparing the 
left and right hand sides shows that the scaling dimensions ^n\p should gather into families with 
^n;p = 'Cn;0+p(2— C) with the modcs ipn;p related by the descent equation ( [7^ ) and dtipn;p = — 
The modes ipn;p are called slow modes because their expectation values, as defined in eg. ([7l|) , 



increase with time slower than what the generic hyper diffusive (72) behavior would have predicted 



Expansion ( [751) may be interpreted as a kind operator product expansion. 

4.4 Breakdown of Lagrangian flows. 

Limit of coincident initial points. 
Ill-definiteness of trajectories. 

Let us now see what happen to -Pk=o initial points are approaching each others. As the 
velocity field is time reversal invariant, this is equivalent to approaching the final points. If the 
Lagrangian trajectories were well-defined, ie. if the Lagrangian trajectories were uniquely defined 
by their initial or final positions, we would have: 

P^^\x,t\y,s) =5ixi-X2)P''^-^\x,t\y,s) (76) 

yi=y2 

Instead for ^ < 2 the limiting distribution, 

pS(x, t\0, 0) = \eM-tM%)] (77) 

is a regular function, ie. it is not concentrated on the delta function. This may be checked in the 
limit ^ — > in which M% becomes the Laplacian operator. These functions give the probability for 
trajectories starting at the same initial point to split. A more precise description of this splitting 



M9, = D'z- 



is provided by the expansion ([75D . 

The distribution of two particles starting at identical initial positions may be computed exactly 
in the rotation invariant sector. After changing variable from r to z = r^^"^)/^, the differential 
operator becomes 

d_ 2 h{h - ly 
'^dz' ^ z2 

with h = ~ \ a-iid D' = D{d — 1). The distribution of two particles P^5q(2;, t|0) is then solution 

of {dt + M'^)pI^% = with the boundary condition Pi%{z,t = 0|0) z^'^dz = dz6{z). The result 
is: 

2 7 

P'ilair, t|0) r'^-^dr = const, {^f exp (-z^/AD'i) (78) 

The violation of equation (^) implies a breakdown of the Lagrangian flow in Kraichnan's 
velocity fields. This breakdown is linked to the irregularities of the velocity; irregularities which 
imply the existence of many solutions to the Lagrange equation (^). It is similar to the one 
predicted by Richardson's law. A meaning to the distribution of these solutions at fixed velocity 
realization has been given in |25]. 



This spreading of the trajectories implies an information loss, as knowing that two particles are 
at identical positions at time to does not guarantee that this information will be true at later time. 
A more quantitative formulation of this information loss will be welcome. 
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4.5 Batchelor limit. 



Trajectory pdf's for smooth fields. 
Trajectories are well-defined. 
Lyapunov exponents and universality. 



The Batchelor limit is the limit of smooth velocity field, ie. ^ = 2 in Kraichnan's model. In 
this limit Lagrangian trajectories have different behaviors, closer to those of trajectories in chaotic 
systems. 

The probability distribution function for two trajectories satisfies {dt + Ai2)P^'^^ = which for 
^ = 2 and in the rotation and translation invariant sector reduces to: 



) P''>(r,(|r„,0) =0 



r'^ ^ dr dr , 

with the boundary condition r'^~^P2{r, 0|ro, 0) dr = 5{r — tq) dr. The solution is: 

P(2)(^,i|^0,0) = -^Ll=exp (^-^ (log(r/ro) - dD'tf^ (79) 

with D' = D(d — 1). Its limit of initial coincident points is such that: 

lim P^'^'^{r,t\ro,0)r'^-^dr = 6{r)dr 

ro^O 

Thus eq.(|7^) is satisfied meaning that trajectories are well-defined. This could have been expected 
since for ^ = 2 Kraichnan's velocity field is smooth because the typical value of its variation between 
points distant by r is 5ti ~ r. 

However the trajectories separate exponentially as it can be seen from eq.(^) or from eq.(72) 
in the limit ^ — > 2. This exponential separation prevent us to use the Batchelor limit to cure the ill- 
definedness of Lagrangian trajectories. Indeed assume that the turbulent velocity fields have been 
regularized such that they become smooth on scales less than say rj but remain irregular above this 
scale. Lagrangian trajectory whose initial points are distant less than rj first enjoy an exponential 
separation until they are distant by a length of order r], and then follow different trajectories since 
they are well separated. 

In this smooth limit, trajectory statistics is solvable by reducing it to a group theoretical 
problem |^, 27]. Let us just illustrate this point by computing the Lyapunov exponents of the 
Lagrangian fiows [p^ j. Since the velocity is smooth, the velocity variation may be parameterized 
as 6u"'{x,t) = a°'^{t)x^ with a{t) a traceless d x d matrix, a{t) G sl{d). For the Kraichnan's model 
a{t) is gaussian with two-point function: 

(a^^t) a'"^{s)) = 2D ((d + 1)6'"'5'"^ - S^'H"'^ - 6'"^6^^) 5{t - s) = C''''"^ 6{t - s) 

The equation for the separation of Lagrangian trajectories is then: 

±{t) = a{t) ■ x{t) (80) 

Its formal solution is 

x'^(t) = G"^(t)x^ , G{t) = P exp (^J^^ dsa{s)^ (81) 

The matrix G{t), which belongs to the group SL(d), detG = 1, codes the information for the 
Lagrangian flows. The later may thus be seen as a process on the group SL{d) with 

G{t) = a{t)G{t) (82) 
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It is a Markov process which is fuhy determined by the probability transition Pt{G,Go) from the 
group element Go to G. It is formally defined as Pt{G,Go) = {6{G^^G{t,to))) with 6{-) the Dirac 



measure for the Haar measure on SL(d) and G(t,to) solution of eq.(82) with initial condition Go- 
The associated Fokker-Planck equation, which may be derived using techniques explained above 
for the Navier-Stokes equation, reads: 

d, - Ig"'''^'^ Lab) Pt{G, Go) = (83) 

with the initial condition i-t(G, Go)|t=o = S{G~^Go). Here Lab are the vector fields corresponding 
to the infinitesimal left-action of gl{d) on d x d matrices, Labf{G) = 4: f {e^^''^ G)\i=Q with Eab the 



elementary matrices. The Fokker-Planck equation (83) may alternatively be written as sum of 
Casimir operators: 



( 



dt - D{d + 1) J2 - DdH'^] Pt{G, Go) = 



with H"^ = J2ab LabLba and J"^ = \ Y^abi^ab - Lbaf where Lab = Lab-\ J2a Laa are generators of 
sl{d). 

Rotations of the initial (final) points correspond respectively to the right (left) actions of SO{d) 
on SL{d). So observables invariant by initial and final rotations are functions on the coset space 
SO{d)\SL{d)/ SO{d). Examples of functions on this space are provided by traces tr(GG^) in any 
representation of SL{d). It is convenient to Iwasawa decompose G as G = OiN with Oi G SO{d) 
and N a lower triangular matrix. The later may be factorized as = DO2 with O2 G SO{d) 
and D diagonal, D = diag(zi, Z2, - • • , z^) with Hj = 1- This identifies points of the coset space 
SO(d)\SL{d)/SO{d) with diagonal matrices. The eigenvalues Zj describe how blobs of initial 
particles are stretched. They typically increase exponentially. So that at large time the logarithms 
of these eigenvalues are asymptotically linear in time: 

log Zj ~ Xjt as t — > 00 

The coefficient Xj are called the Lyapunov exponents. 

On SO{d)\SL{d)/SO{d), the Fokker-Planck equation iM) reduces to: 



dt - DdJ2dl + DiJ2 dyf - DdY^cothivi - yj)dy^ Pt{y; yo) = (84) 
j j i+i j 

with Zj = expyj, J2jyj = 0- This is the Schroedinger equation for a Calogero-Sutherland type 
hamiltonian. It is also the Fokker-Planck equation for the Langevin dynamics: 

dtyj = DdY^ coth(yi - yj) + (85) 

where are random gaussian variables with zero mean and {£,j{t)(,i{s)) = 2D{d5ij — l)(5(t — s). 
The typical dynamics is simple to analyze. Suppose that at initial time the eigenvalues are ordered 



Ai > A2 > • • • > Arf. The interacting term in (|85D preserves this order and increases the splitting 



between the eigenvalues. So at large time Ai » A2 3> • • • 3> A^ and we can then approximate 
coth(yj — yj) by ±1. The effective dynamics of the eigenvalues then reads: 

dtyj=Dd{d-2j + l)+ij 

This shows that the Lyapunov exponents are \j = Dd{d — 2j + 1). This effective dynamics also 
shows that at large time the statistics of the logarithm of the eigenvalues is gaussian. 
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Trajectories in smooth but time-correlated velocity fields have been analyzed in [29|. The 
Lagrange equations are then eq.(|80|), or equivalently eq.(|^), but with a{t) time-correlated. De- 
composing the group element as above G = ON with O G SO{d) and triangular, the time 
evolution reads 0~^0 + NN^'^ = 0^^a{t)0. Let as above Zj = expy^ be the eigenvalues of N, 
then ZjZ~^ = {0'^aO)jj = ajj. This is simply integrated as: 



yj{t) = yoj + [ dsajj{s) 
Jo 



If r is the correlation time of a, the above integral is a sum of t/r independent variables. We 
can apply the central limit theorem for t large. The bulk of the distribution of the y'jS is thus 
gaussian with mean average increasing linearly with time. The proportionality coefficient \j are 
the Lyapunov exponents with 

1 r 

The gaussianity of the statistics illustrates the fact that the statistics of Lagrangian trajectories 
in time-correlated but smooth velocity fields are described by effective universal dynamics. See 



ref.[29| for a discussion of consequences of this property on scalar turbulence. 



4.6 Generalized Lagrangian flows and trajectory bundles. 

Non-uniqueness of trajectories and parameters as random variables. 
Generalized Lagrangian flows. 

Prom the two preceding sections we learn that Lagrangian trajectories are well-defined for ^ = 2, 
ie. when the velocity field is regular, but ill-defined for < 2 in the sense that eq.(^) is broken. This 
arises from the fact that there is many solutions to the equation x = u{x, t) if u{x, t) is not smooth 
enough. To exemplify this fact consider x = x'^ with < cr < 1 which possesses two solutions with 
initial conditions x(0) = 0, namely x{t) oc t^/^^'^ or x{t) = 0. So the probability distribution p(^) 
are actually not describing discrete set of particles but diffuse clouds of trajectories. 

More generally, let us parameterize the set of solutions to Lagrange equation as x{uj, t) so that 

x{ui,t) = u{x{ui,t),t) (86) 

with u{y, t) some realization of the velocity field. One may thought of uj as labeling the initial 
condition plus the extra parameters needed to specify the solutions. To make it more precise, but 



more formal, one assumes [30, |3l| that uj ft with Q a space of events equipped with a probability 
law, ie. equipped with a measure fi{du!). The data of map x : Q x [ti,tf] —>■ whose output is 
x{ui,t) has been called generalized Lagrangian fiow. 

Incompressibility may be forced by demanding that the measure fi{duj) satisfies: 

fi{du;)ipt{x{uj,t))= dyiptiy) (87) 



for sufficiently regular functions (pt{y) on R . Indeed time differentiating eq.(|87D leads to the weak 
incompressibility condition / dyu{y,t) • Vipt{y) = 0. 

The case of well-defined trajectories corresponds to ^2 = R"^ with d the dimension of the ambient 
space and fj, the Lebesgue measure. The parameter oj is then the initial position. In the opposite 
case with a breakdown of the Lagrangian flow we should talk not about single trajectory but only 
about what may be called a trajectory bundle with some probability distribution law, see flgure. 
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For example the probability to find a trajectory in a domain D C R'^ a time t knowing that it 
was inside a domain T>q C R'^ a time to is given by: 



Pt^toi'D ,'Do\u) = dx ii{duj)S{x - x{uj,t)) 



wih J^Do,to = {lj G s.t x{ijo,tQ) £ T>o} the subset of events with trajectories initially inside T>q. It 
is worth stressing that this is a probability distribution at fixed velocity realization. 

In the case of Kraichnan's model this probability distribution for a single trajectory bundle is 



Ptfy,{'D,'Do\u) = dx dy R^=o{x,t\y,to) 

JV JVo 



with RK=Q{x,t\y,tQ) = P • exp /j(,(V-n)(s) defined in eq.(41,49). This is still a random func- 
tion as Kraichnan's velocity is itself random. The average of it gives the probability distribution 
Pi=;'\x,t\y,to). 

It will be quite interesting to extend the previous analysis (conserved and slow modes, ill- 
definedness and probability distribution functions of trajectories, universality, etc..) to non-smooth 
and time-correlated velocity fields. 



32 



5 Burgers turbulence. 



Burgers equation and shocks. 
Bifractality and velocity pdf's. 



The forced Burgers equation in 1 + 1 dimension is the following: 



dtu + udxU — ud'^u = f 



(89) 



The force is supposed to be gaussian with covariance {f{x,t)f{y,s)) = C{x — y)5{t — s). There is no 
pressure and thus no incompressibility condition is imposed on the velocity field u{x,t). Although 
its properties are very different from those of the Navier-Stokes equation both equations possess the 
same non-linearity. The Burgers equation hence provides an interesting toy model for turbulence 
in which progress have recently been done, cf eg ^] and references therein. The main point 
about Burgers turbulence is that one knows the geometrical structure responsible for intermittency: 
these are shocks. 

Solutions of eq.(p^ developp shocks in the inviscid limit. The occurence of these shocks are 
esay to understand by considering the Burgers equation in absence of force. In the inviscid limit it 
reduces to the Euler equation dtu+udxU = whose solutions are such that u{x, t) = uo{x — tu{x, t)). 
Assume for example that the initial condition ^0(2;) is such x = —auQ{x) — tuo{x) with a and r 
positive. Then at latter time u{x,t) satisfies x = —au^{x,t) — (r — t)u{x,t). Hence u{x,t) would 
be multivalued for t > t which means that a shock has developped at time r. 

The Burgers equation may be rewritten in different ways. First setting u = dxh and / = —dxV 
it becomes the so-called KPZ equation : dth + \{dh)'^ — vd^h = —V. Furtheremore, setting 
h = — ^log'i/', it becomes equivalent to the one dimensional Schrodinger in a random potential: 



For typical realization the velocity profile will be piecewise smooth, meaning that the profile is 
made of a succession of intervals, on which the velocity is smooth, these intervals being separated 
by shocks. 



More explicitly as a function of the position the velocity is discontinuous so that its derivative 
may be decomposed as: 




u(x,t) 




Figure 6 : Typical profile in Burger's trubulence. 




(90) 



s 
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The sum is over the shocks, which are known to form a discrete set. The points Xg are the shock 
locations and their amphtudes. The function ^{x, t) is regular except at the points Xg where it 
may have discontinuities. So dx^{x,t) may also have delta function singularities. 

The positions and the amplitudes of the shocks are time dependent. The velocity of the shocks 
is the mean value of the velocity on the two side of the shock so that 

Xg = Vs with u{xf ,t) = Vs ± as/2 (91) 

By definition the amplitude of the shock is Ug = u{xf) — u{x^). The shock dynamics follows directly 
from Burgers equation. Time derivating eq.(pl|) yields dt{vs ± ^) = Xg{dxu){xf) + {dtu){xf ). 
But outside shocks there are no anomaly in the inviscid limit and the inviscid Burgers equation 
dtu + udxU = f holds. Applying it on the two sides of the shocks, ie. for x xf, gives {dtu){xf ) + 
(vg ± ^){dxu){xf) = f{xg). By definition {dxu){xf) = £,{xf ) so that: 

Vs = -^{axt)-a^j)) + f{Xg) 

as = -y(C(4)+e(x7)) 

Time evolution of higher derivatives of the velocity at a shock may be derived in a similar way. 
Namely, let = ^{d^u){xf) and /„ = ^{d^ f ){xg), then applying Burgers equation on the two 
sides of the shock gives: 

i^-{n+l)vgCt+i+ E P^P^^ = fn 
p+q=n+l 

At finite viscocity there is no shock, ie. no discontinuity of the velocity. The shocks are becoming 
domains of width of order l/v over which the velocity decreases rapidly. Inside the shock domain 
located at Xg the velocity will be of the form 

Ug{x,t) = w{{x — Xg)/i>,t; i^) with 'w{z,t;i')='Wo{z,t) + i"Wi{z,t) + --- (92) 

The condition to connect this solution to the rest of the velocity profile on the two sides of the 
shock is that w{±oo) = u{x^). Solving the Burgers equation with this ansatz gives to lowest order 

wq{z) = Vg — ^tanh(asz/2) (93) 

This knowledge may for example be used to illustrate the presence of dissipative anomalies. 
Consider the dissipative field e{x) = ^{dxu)'^. For a given realization, u{x) is smooth outside 
shocks and u{dxu)^ vanishes there. Thus only shocks contribute to the dissipation. To find the 
contribution of a given shock, localized at Xg, one uses the representation (^) of the velocity inside 
the shock to compute the dissipation. One has ^{dxUg)'^ ~ ^-^^6{x — Xg) as ^ 0, so that e{x) is 
the sum of the shock contributions: 

<^)= E ^^^x-xg) (94) 

sishock 

More generally, one may also consider family of operators defined as products of the dissipation by 
functional of the velocity, eg. e(x) 5{v — n(x)), or consider higher derivatives of the velocity. They 
are all localized on shocks and they code for the structure of the velocity profile inside and around 
the shocks. 
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Shocks dominate many properties of Burgers turbulence. The first consequence is that the veloc- 
ity structure functions are bifractal so that Burgers turbulence manifests a pronounced intermittent 
behavior. One has for x small and positive: 

with ps the shock density. In other words, for < p < 1 the moments are dominated by the smooth 
parts of the velocity, whereas for 1 < p they are dominated by the shocks. This may be understood 
as follows. If there is no shock in the interval [y,y + x], for small x the variation of the velocity 
between the two points y and y + x is 6u{x) ~ with ^ the derivative of u. The probability for 
such an event is (1 — psx) for small x. If there is a shock in the interval [y,y + x], the velocity 
variation is 5u{x) ~ with the amplitude of the shocks. The probability for this event is psX. 
Thus the moments may be evaluated as: 

{{5unx))^{i-p,x){{^xr) + iPsx) K) 

It is the first term which dominates for < p < 1 and the second for 1 < p. This reproduces 
eq.®. 

Another interesting quantity is the probability distribution function for the gradiant of the 
velocity. At finite viscocity it is formally defined as P{C) = {S{C — dxu)). Its inviscid limit is the 
probability distribution function of £,{x,t) the regular part of dxu{x,t) as defined above. As usual, 
the stationarity condition leads to a generalized Fokker-Planck equation: 

CP(C) + d^iePiO) + Ci a|P(C) = ^^((a^n) 6'{C - dxu)) (96) 

with Ci = — ^C"(0). The r.h.s. is the analogue of the dissipative anomaly but for higher gradiant 
of the velocity. It is no vanishing due to the presence of shocks and only shocks contribute to it. 

Universality appears in the tails of this probability distribution function in the inviscid limit. 
Consider first the tail for large positive vecolity gradiant, Q +oo. It is the regular part of 
the velocity which dominates in this domain because the velocity gradiant is large and negative 
inside shocks. So for C, +oo it is thus legitimate to neglect the anomalous l.h.s. in the Fokker- 
Planck equation (p^). The differential equation may then be solved and, with appropriate boundary 
conditions, it gives |Q: 

P(C) =i const. Ce-'^'/^'^i for C ^ +oo (97) 

The behavior of -P(C) as C — > — cxd is more tricky because shocks contribute significantly in that 
domain. One expects a power law behavior: 

P(C) ^ const. |Cr° for C ^ -oo (98) 



A nice analysis presented in shows that realizability constraints impose a < 3 while a finer 
analysis of shock formations yields to a = 7/2. 

It will be interesting to have a statistical description of the forced Burgers turbulence as precise 
as one has for the decaying Burgers turbulence, see eg. [32, 33]. 
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